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When an ion confined in a linear ion trap interacts with a coherent laser 
field, the internal degrees of freedom, related to the electron transitions, couple 
to the vibrational degree of freedom of the ion. As a result of this interaction, 
quantum dynamics of the vibrational degree of freedom becomes complicated, 
and in some ranges of parameters even chaotic. We analyze the vibrational 
ion dynamics using a formal analogy with the solid-state problem of electron 
localization. In particular, we show how the resonant approximation used in 
analysis of the ion dynamics, leads to a transition from a two-dimensional (2D) 
to a one-dimensional problem (ID) of electron localization. The localization 
length in the solid-state problem is estimated in cases of weak and strong 
interaction between the cites of the 2D cell by using the methods of resonance 
perturbation theory, common in analysis of ID time-dependent dynamical 
systems. 
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INTRODUCTION 



The problem of quantum dynamics for Hamiltonian systems with time-periodic pertur- 
bation (TPS) can be formulated in terms of an equivalent solid-state problem of electron 
localization on a lattice. Such kinds of connections were discussed for different models in&i 
(see also references therein). However, most results are obtained for quantum kicked systems, 
such as a quantum kicked rotor or a quantum kicked oscillator. These systems are convenient 
for both analytical and numerical analysis because instead of differential equations one can 
use discrete quantum maps. At the same time, the existence of periodic kicks suggests that 
the external field involves an infinite number of harmonics with equal amplitudes. In more 
common physical situations, there are only few harmonics in the perturbation. In particular, 
such a situation occurs when an ion trapped in a linear ion trap interacts with two laser 
field with close frequencies!. In this case, the internal degree of freedom of the ion (related 
to the electron dynamics) interacts with the vibrational degree of freedom. This interaction 
can result in complicated and even chaotic dynamics of the vibrational degree of freedom 
of the ion. The analysis of the stability of the ion in this system can be performed using a 
model of a quantum harmonic oscillator perturbed by a monochromatic wave (MPO)i. 

In this paper we show that the problem of stability of a MPO can be formulated in terms 
of localization of an electron in a 2D solid-state system (SSS). The resonance approximation, 
common in treatment of TPSs, is used to reduce the effective dimensionality of the SSS in 
the case of relatively small interaction of the ion with the laser field. In order to compare two 
completely different systems, a similarity in the formal description of the TPS and a space- 
periodic SSS is exploited. Namely, in the TPS we use the time-periodicity of the perturbation 
and employ a Floquet formalism, while in the SSS we exploit a space-periodicity and use 
the Bloch theorem. 

The paper is organized as follows. In Section I we describe the MPO model of an ion 
trapped in a linear ion trap and interacting with two laser fields with close frequencies. In 
Section II, we discuss the general procedure which allows to treat an ID TPS on the same 
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ground as a 2D SSS. In the case of a small perturbation, the resonance approximation is 
used in Section III to decrease the effective dimensionality of the SSS from two to one. 
The localization length in the SSS is estimated in Section IV by calculating the size of the 
chaotic region in the corresponding TPS, in the situation when the interaction between the 
sites of 2D SSS is strong. Concluding remarks are given in Conclusion. 



I. THE VIBRATIONAL HAMILTONIAN 

In the following!, we assume that two laser beams, designated the pump (p) and the 
Stokes (s), with slightly different frequencies, u p and u s , respectively, interact with the ion 
trapped in a linear ion trap. Both beams are assumed to be plane polarized in the z-direction 
with the amplitudes of the electric field, and £^ s \ and the wave vectors, k p and k s . The 
Hamiltonian, including the effect of the harmonic evolution of the ion along the weak axis 
of the trap (but excluding the internal free evolution) is, 

where where p and x are the x-components of the momentum and the coordinate of the ion, 
t is the time, M is the mass of the ion, u is the frequency of the ion vibrations in the linear 
trap, e = 2\k Sf^'Ef'* , Cl = cu p — cu s , k = (k p — k s ) • e x , e x is a unit vector in x-direction, 
X = A7reo/4i/ 3 A [v and A being, respectively, the wavenumber and the Einstein A coefficient 
for the transition between the upper and lower manifolds, A is the laser detuning and eo the 
permitivity of free space). 

In the dimensionless form the Hamiltonian ([!]) reads, 

" = = 4 + T + f C08( - Y " " T) = + V{X ' T) ' (2) 

where Ho is the Hamiltonian of a linear oscillator, 

X = kx, T = ut, e =—— , h = ——, fi=- = N + 6. (3) 

Muj z Mu uj 
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Here h is a dimensionless Planck constant, N is the (positive integer) resonance number, 
and 5 is the detuning from the resonance. 

The classical analog of the Hamiltonian (|2]) is, 

X 2 P 2 

H = — + — + ecos(X-fiT), (4) 

where P = kp/Mu is the dimensionless momentum. 

In the action-angle variables, (/,#), the classical Hamiltonian (0) takes the form, 

H = I + ecos[kr(I) sintf — fir], (5) 



where X = kr(I) sin$, P = kr(I) cos$, kr = y/X 2 + P 2 = y/2I is the dimensionless 
amplitude of oscillations, I is the dimensionless action, measured in units of Iq = Mu/k 2 , 
and d is the phase of oscillations. 

II. CONNECTION WITH A 2D SOLID-STATE LOCALIZATION PROBLEM 

We write the solution to the Schrodinger equation, 

it^l = ^ { X,r), (6) 

in the form of series over the eigenfunctions, \n) = n (X), of the harmonic oscillator Hamil- 
tonian, Hq, 

oo 

9(X,T) = 52c n {T)\n). (7) 

n=0 

Then we obtain the equations for the complex amplitudes, c n (r), 

dc m (r) 00 
ih — ^ — = h(m + l/2)c m (r) + e ^ (m\ cos(X — //r)|m + n)c m+n (T) — (8) 

^ T n=-m 

h(m + l/2)c m (r) + l £ (e-^F m , m+n + e^F^ m+n )c m+n (r). 

n=—m \ / 

In Eq. (H), -F m m+n is the matrix elementi, 

i n h n/2 e~ h/A ( h\ 

F m , m+n = (m\e iX \m + n) = L n m - , (9) 

2"/ 2 J(m+l)(w+2)...(m + n) V 2 / 



where is the Laguerre polynomial. When m 3> 1, the Laguerre polynomials can be 
expressed in terms of the Bessel functions, J n i, as 

^(^(x)^^' (10) 
where the argument of the Bessel function, y2mh = kr m , is the quantized dimensionless 
amplitude of oscillations of the harmonic oscillator. Using Eqs. (P) and fliTf ) the matrix 
elements can be written in the form, 

F m ,m+n = i J n (V2mh). (11) 

J(m + 1) . . . (m + n) 

Since the Hamiltonian (|2]) is periodic in time the solution of the Schrodinger equation 
(|) can be written as, 

cl(r) = e-^AKr), (12) 

where a q is a quasienergy, measured in units of Mu 2 /k 2 , c-^(r) is the quasienergy (QE) 
function, and A^t) is a periodic function with the period, T = 27r//x, 

Al(r + 2n/^ = Al(r). (13) 

The quasienergy functions are the eigenfunctions of the evolution operator, U, for one period, 
T, of the external field, 

U(T)cI(t) = e^ h cl(r). (14) 

The evolution operator for one period of the external field, U{T) is, 

U(T) =fe- i ^ {l{T)dT , (15) 

where T is the ordering operator, H(t) in our problem is given by Eq. (|2|). In our numerical 
calculations, presented below, we consider only the QE states at the time r = 0, so that 
4,(0) = Al(0) = Alft 



1 If the spectrum, a q , and QE functions, A^, are known, one can trace evolution of the quantum 
system at the times t s = sT, where s = 0, 1, 2 . . . (see, for examplei). 




Substitution flT2|) in (§) gives the equations for A q m , 

<V^(vO = ~ %h ^0^ + h ( m + V2)^)+ (16) 
£ g ( e -^F m>w+n + e^F* )m+n )^ +n (^), <p = (jlt. 

n=— m \ / 

Expanding the function A^tp) in a Fourier series, 

oo 

KSp) = E (17) 

we derive the following equation for the amplitudes, A q ml1 



e 00 



EqAm,l — h( m A^)^m,J + 2 E I ^m,m+n^m+n,l+l + ^m,m+n^m+n,l-l )' (18) 



n=— m 



where .E^ = a q — h/2. Equations (|18|) can be interpreted as a problem for electron localization 
on a 2D lattice. Indeed, one can consider the complex coefficients, A q mh as the complex 
amplitudes of probability to find an electron on a 2D lattice at the site (m, I), where < 
m < 00, — 00 < Z < 00. Some possible transitions in the system (0) are shown in Fig. 1 (a). 
The particle, initially located at the site with the indexes (m , lo) can jump at the sites 
(mo ±n,/oil), where n is an integer number. 



III. THE RESONANCE APPROXIMATION 



When the interaction amplitude is small, e < 1, the 2D SSS described by Eq. (|T8| ) can 
be reduced to the ID SSS described by the equation with only one index. We divide both 
parts of Eq. (|18|) on /i and, taking into account that 5 <C N and ~ 1/N — 5/N 2 , obtain 



E n .„ , (m 5m ,\ .„ e x 



We assume that 5 is small, so that 5m/ N 2 <C 1 for all considered values of m, or 5 = 0. 
Then, in the zeroth order approximation we have from Eq. (O) , 



m A .„ 



= -fr^nj- (2°) 



It follows from Eq. (0) that if A q ml ^ 0, than E^f/ih = (m/N) - I. Since the ratio, 
EfP> I fj,h, is defined by modulus 1 (see Eq. (|1|)), we can write, Efj[ih = 0Q Then Eq. (H) 
takes the form, 



(m - Nl)A q = 0. (21) 



Hence, 



A q m l = for m ^ iVZ, (22) 



<,=^ for m = Nl. (23) 
The next order approximation for m = Nl yields, 

(E q - h5m/N)Al = ~{F m , m+N Al +N + F^ m _ N A q m _ N ), (24) 



2 We assume that (m/N) — I is an integer for some initial state tuq. If not, one may intro- 
duce the quasienergy, E' = E q — [ih{mo/N} and a new index, m! = m — {mo/N}N , and solve 
Eqs. (|i~9| ) and (^) for ^ and m! instead of E q and m. Here {x} is the fractional part of x. 
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where E q = E^> (we do not consider the higher order approximations). As one can see 
from Eq. (]2l|), in the m-direction only hops on the distance N are allowed. Thus, the 2D 
problem, given by Eq. ([18]), is reduced in the case e C 1 to the ID problem described by 
Eq. (H. 

The localization properties of the quantum states in the resonance approximation, given 
by Eq. (|24|), are defined by the structure of the matrix elements, -F mjm +Ar. If the matrix 
elements are periodic functions of m, all the eigenstates are extended and the spectrum 
is continuous. This situation is common for solid-state systemsi. In the system under 
consideration, the matrix elements (|TT| ) are non-periodic. On this reason, as will be shown 
below (see also RefsBi) the quantum states are localized and the spectrum is discrete. 

The matrix elements given by Eq. (|IT1) , oscillate as a function of m. At the points, m , 
where the matrix elements are close to zero, 

J N (v / wT) « 0, (25) 

the transition probability is very small. As a consequence, such points becomes the dynam- 
ical barriers to the probability and divide the Hilbert space, labeled by index, m, 
into the relatively independent parts — resonance cellsl. Most of the eigenstates given by 
Eq. ([24]) are concentrated inside these cells. The average localization length, (Aj), for the 
states in the i-th cell does not exceed the size of the cell. The cell boundaries are defined 
by Eq. (|25|), i.e. (A,) < m i+ i — rrii, where and m i+1 satisfy Eq. (p5|) , so that \J1mih and 
\f2m i+ ih are, respectively, the i-th and i + 1-th roots of the Bessel function in (p5|). 

Some characteristic QE functions given by Eq. (p4|) are illustrated in Figs. 2 (b) - (e), 
for small e. The boundaries of the resonance cells are marked by arrows. One can see from 
Figs. 2 (b) - (e), that the eigenfunctions are localized inside the cells, but, on the other 
hand, each eigenfunction is delocalized over m inside a single cell. For example, for the 
initial states in Fig. 1 (b) with m = 72 the transitions will occur in the region 66 < in < 176 
(inside the first cell in Fig. 2 (a)). Note, that for small values of e and when 5 = 0, the 
localization properties of our system are independent of e. This means that an arbitrary small 
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perturbation, e, initiates transitions between the sites on the effective rectangular lattice. 
For small values of e, these transitions take place on ID sub-lattices shown in Fig. 1 (b). 
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FIG. 2. (a) The matrix elements (in arbitrary units), F m)7n+ 2, and (b)-(e) some characteristic 
QE functions given by Eq. (|24| ) with h = 0.2, e = 0.02, 5 = 0, N = 2. Only even values of m are 
included. 



Except for the localized (in the resonance cells) eigenf unctions there exists few eigen- 
functions which are delocalized over several resonant cells (DF). One of these representative 
eigenfunctions is shown in Fig. 3. The DFs have maxima in the regions near the boundaries 
of the cells marked in Fig. 3 by arrows. Thus, if the initial state is located near the boundary 
of a resonance cell, say, at m = 66 in Fig. 3 (in the region near the first arrow), then this 
state will propagate at a large distance in m, over the ID sub-lattice. This distance can be 
much larger than the size of the single resonance celfi 
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FIG. 3. The characteristic DF (at even m),q = 323, N = 2, h = 0.2, e = 0.02, 5 = 0. 
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FIG. 4. The plot m, versus A q for iV = 2, h = 0.2, e = 0.02 and (a) <5 = 0, (b) <5 = 0.001. The 
boundaries of the resonance cells are marked by the arrows on the m 9 -axis. 

The structure of the eigenfunctions can be better understood from the plot of the mean, 
m q — J2 m l^ml 2m > versus variance, A q = [J2 m \^m\ 2 ( m ~ m q) 2 ] 1 ^ presented in Fig. 4 (a). 
Each eigenf unction, A q m , is represented by one point in the figure. One can see that most 
of the eigenfunctions are localized inside the resonance cells since their means are located 
inside the cells and their variances do not exceed the size of the cell. Each row on the figure 
is formed by the eigenfunctions of one cell. If the initial state is located inside the resonance 
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cell, the eigenfunctions of this particular cell define the quantum dynamics. These states 
make the quantum dynamics localized inside the cell and, at the same time, delocalized over 
the states, m, inside the cell. In the corresponding solid-state model (|24]), the localization 
length at small e can be identified with the size of the resonance cell. 

The DFs are represented in Fig. 4 (a) by the scattered points with large variances. One 
of the DFs, marked in Fig. 4 (a) by an arrow, is shown in Fig. 3. The DFs can not be 
attributed to a definite resonance cell since their variances are larger than the size of single 
cells. As a consequence, the DFs cause derealization of the states initially concentrated 
near the boundaries of the cells. However, as shown in Ref.@, the localization length remain 
finite, because the matrix elements ( pTT|) are non-periodic and their amplitudes decrease with 
m increasing (as m" 1 / 4 at m > 1). 

In the case when the detuning from the resonance (see (||)) is not equal to zero (8^0), 
the character of localization depends on the position of an initial state, mo (see Fig. 4 (b)). 
In the region m 3> m max = eN / h5 all the states remain exponentially localized in m, since 
in this case Eq. (|24|) has the solution, 

E q = (h5m/N)5 mtq , A q m = S m , q . (26) 

If mo <C ni max the above discussed effect of localization over the resonance cells takes 
place. In the intermediate case, when m > m max , the character of localization depends 
on the position of the state m inside the resonance cell. (For the parameters in Fig. 4 (b) 
m max — 200.) If m is located near the boundary of the resonance cell where the condition 
(p5p is satisfied, Eq. (|4]) has the localized solution (p6"l). Most delocalized functions have 



their mean, m q , at the center of a resonance cell. 

As follows from Fig. 4 (b), in the region m > m max the DFs are absent, since the 
variance of each function is much less than the size of the cell, whose boundaries are marked 
in Fig. 4 (b) by arrows. Moreover, the variances of the eigenstates in the near-resonance 
case in Fig. 4 (b) are substantially smaller than the variances in the exact resonance case 
shown in Fig. 4 (a). Hence, at small e, the increase of the value of the detuning, 8, always 

11 



leads to localization of the quantum states in the discussed model. 

Most of the localization properties of the eigenfunctions, given by Eq. (|24]), are the 
quantum manifestation of the classical behavior in the phase space. The classical phase 
space in the variables (kr(I), 6)), where 6 = N$, mod 2n, generated by the exact classical 
Hamiltonian (|5|), is shown in Fig. 5 (a) for the exact resonance case (5 = 0) and in Fig. 5 (b) 
for the near resonance case (5 = 0.001). 
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FIG. 5. The classical phase space for (a) 5 = 0, (b) 5 = 0.001; N = 2, e = 0.02. The boundaries 
of the resonance cells are marked by arrows. 

As one can see from Fig. 5 (a), in the case 5 = the classical phase space is divided 
into the resonance cells. (Fig. 5 (a) shows only the first seven cells.) The boundaries of 
the cells, kr(Ii) = y/2Ti, are marked in Figs. 5 (a), (b) by arrows. As shown in Refs.ii, 
in the quasiclassical limit the i-th boundary of the classical cell, kri, in Figs. 5 (a), (b) 
corresponds to z-th boundary of the quantum cell, m;, in m q - axis in Figs. 4 (a), (b), so 
that kr~ ■■ 



s /2T i = y/2hm~i. 

Each row of points in Fig. 4 (a) is formed by the eigenstates responsible for the dynamics 
in the corresponding quantum cell. From comparison with the classical dynamics in phase 
space we can now describe the localization properties of the quantum states, discussed 
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above. Each value of m in the quantum system corresponds to a quantized classical action, 
J m = mh, or to the quantized dimensionless oscillation amplitude, kr m = \j2mh. Each value 
of action I m (or kr m ) corresponds to the set of classical trajectories. Moving along some 
classical trajectory the particle with some initial value of action, J m , can accept the other 
values in the interval I mi < I m < I m2 . The corresponding eigenstate will be delocalized 
over the unperturbed states with the numbers m in the interval, m\ < m < ni2. From 
the form of the trajectories in Fig. 5 (a) one can see that in the case of exact resonance 
all quantum states of the single quantum cell should be delocalized over the resonance cell, 
since in the phase space both the extremal values, I mi and I m2 , which limit the resonance 
cell, can belong to the same trajectory. 

Similar arguments can be used to analyze the quantum-classical correspondence in the 
near resonance case. As follows from Fig. 5 (b), at 5 ^ in the phase space there is only 
a finite number of resonance cells (two cells in Fig. 5 (b)). Thus, there is a finite number 
of quantum resonance cells in the Hilbert space in Fig. 4 (b) (the first two cells). In the 
off-resonant region (3th- 7th cells in Fig. 4 (b)) the degree of derealization of eigenstates 
depends on the position of the state in the cell destroyed by the finite detuning, 5. In 
Fig. 5 (b) the least curved trajectories are located near the separatrices, while the the most 
curved trajectories are located near the centers of the destroyed cells. As a consequence, 
in the quantum model the eigenfunctions in Fig. 5 (b) have smallest variance in the region 
near the separatrices and largest variance near the centers of the destroyed cells. 



IV. THE LOCALIZATION LENGTH IN THE CASE OF STRONG 

INTERACTION 

In the previous sections we considered the dynamics only at small perturbation am- 
plitude, e < 1. At large values of e the dynamical chaos appears in the classical^ and 
quantum@0 MPO. As will be shown below, all quantum states in the chaotic area are de- 
localized over the whole chaotic region. The chaotic dynamics in MPO corresponds to hops 
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in different directions in the SSS, as shown in Fig. 1. By estimating the size of the chaotic 
motion in the MPO we will estimate below the localization length in the SSS. 
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FIG. 6. The classical phase space in the case e = 3. Other parameters are: 5 = 0, N = 2. The 
boundaries of the resonance cells are marked by arrows. 
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FIG. 7. The time-averaged probability distribution, (a) e = 0.02, the averaging has been per- 
formed over 100 realizations in the time-interval r = 5000 — 105000 (only the probability at even 
m is shown); (b) e = 3, where the averaging has been performed over 100 realizations in the 
time-interval r = 500 — 10500. Other parameters are: h = 0.2, 5 = 0, N = 2. The boundaries of 
the resonance cells are marked by arrows. 
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In Fig. 6 the classical phase space is shown for e = 3. One can see that in the first two 
cells the motion is mainly chaotic while in other cells the motion remains mainly regular. 
The quantum probability distribution is illustrated in Fig. 7 (a) for the case of small e and 
in Fig. 7 (b) for the case e = 3. The initial state was taken in the form c m (0) = 5 m , mo with 
m = 30 (in the center of the first cell in Figs. 7 (a), (b)). As follows from Fig. 7 (a), the 
quantum particle can tunnel (see also Ref.&) from the initial (first) cell to other resonance 
cells unlike the classical case, where practically all the trajectories in the phase space are 
confined inside the resonant cells (see Fig. 5 (a)). 



1400 



900 



400 



0.15 



(a) 



(b) 



0.05 



-0.05 



100 



200 A 300 
Aq 



400 



-0.15 



400 



m 



900 



1400 



FIG. 8. (a) The plot m q versus a q and (b) the characteristic QE eigenfunction located in the 
chaotic region for e = 3 and h = 0.2, 5 = 0, N = 2. The boundaries of the resonance cells are 
marked by arrows. 

When e increases, the probability distribution in Fig. 7 (b) increases in most of the 
quantum cells, which corresponds to chaotization of motion in the classical phase space. In 
Fig. 8 (a) we show the plot m q (a q ) for the case e = 3, and in Fig. 8 (b) the characteristic 
QE function located in the chaotic area is illustrated. As one can see from Fig. 8 (a), almost 
all QE states in the area of the first two cells are delocalized over both the cells. In other 
words, the QE states are localized inside the chaotic area (first two cells), but not inside 
the single cells, as in the case of small e in Fig. 4 (a). In the chaotic regime, one can find a 
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quantum particle with equal probability in any unperturbed state, m, inside the chaotic sea, 
independently of the position and the form of an initial state located in this region (see first 
two cells in Fig. 7 (b) and 8 (a), (b)). When e increases, more classical and quantum cells 
become chaotic. This results in increasing the area of derealization of the quantum chaotic 
states. Thus, in the regime of chaos the localization length in the SSS may be identified 
with the size of the chaotic area in the TPS. 



V. CONCLUSION 

The resonance approximation, given by Eqs. (|2l| ) - ([24|) can be interpreted in the following 
way. We can re-write the Hamiltonian matrix in Eq. in the form, 

H m ,l;m+n,V = ^m,J;m+n,J'^.0^,l' + Fn,n+m$l' ,1+1 + ^n,n+m^l' ,1-1- ( 2 7) 

The unperturbed motion, given by the zeroth order Hamiltonian, H^ ll = h(n — fil), takes 
place on the infinite ID energy surface (line) in the 2D space. The character of the motion on 
the energy surface depends on the form of the matrix elements of the interaction potential, 
(to, l\V(X, r)\m', I') = F m j. m /ii. If V(X, r), is linear in X, for example, V(X, r) = eX cos fir, 
and fj, = N, where N = 1, 2, . . ., then the motion is unlimited0 , El. When the interaction 
potential is nonlinear in X, V(X, r) = ecos(X — /rr) at \i — N, or when the period of the 
interaction potential, /i, is not equal to iVx (a distance between the sites of the 2D cell), 
i.e. when 5^0, then, as demonstrated in this paper (see also Ref.@), the quantum states 
are localized. 

It is necessary to note that, as was shown irM the quantum and classical dynamics in 
the chaotic regime is practically independent of the detuning, 5, when e ^> 8. So, the results 
concerning the chaotic dynamics also remain valid in the near-resonance case. 

In summary, the regular and chaotic classical and quantum dynamical regimes are an- 
alyzed in the system of an ion trapped in a linear ion trap and interacting with two laser 
field with close frequencies. This system is modeled using a quantum oscillator perturbed 
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by a monochromatic wave (MPO). It is shown that the problem of dynamical stability in 
this system corresponds to the problem of electron localization in a 2D solid-state system 
(SSS). The resonance approximation is used to decrease the effective dimensionality of the 
corresponding solid-state system. This can be done in the case of relatively small inter- 
action between the trapped ion and the laser fields. In the SSS this case corresponds to 
weak interaction between the sites of the 2D lattice. Increasing the interaction amplitude 
results in derealization of the quantum states over the sites of the 2D cell. The area of the 
derealization in the SSS at strong interaction in the chaotic area may be identified with 
the size of the chaotic sea in the MPO. Our results provide understanding of the mechanism 
of stability of an ion trapped in a linear ion trap. They also allow one to estimate the 
characteristic dynamical regimes of the trapped ion and to choose parameters required for 
dynamical stability. 
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